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INTRODUCTION 

In  response  to  various  inquiries  and  suggestions 
concerning  the  report  "THEORETICAL  STUDIES  ON  THE  PLOW 
THROUGH  NOZZLES  AND  RELATED  PROBLEMS"  (AMP  Report  82. IR), 
sisking  for  supplementary  information,  the  following  report 
was  prepared  ^W  Dr.  K.  0.  Friedrichs  of  the  Applied  Math- 
ematics Group,  New  York  University,  with  the  assistance 
of  Dr.  C  DePrima,  Miss  M.  Johnson,  and  Mrs.  B.  Grossman. 
The  present  supplementary  report  describes  in  detail  the 
construction  of  a  perfect  nozzle  and  contains  some  additional 
remarks  concerning  the  accuracy  and  the  theoretical  back- 
ground of  the  construction  described  in  the  previous  report. 

R.  Co\irant 

Technical  Representative 
Contract  OEMsr-945 
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ON  SUPERSONIC  COMPRESSORS  AND  NOZZLES 

In  Appendix  II  of  the  previous  report  a  general 
method  was  described  to  determine  contours  for  a  con5)ressor 
duct  (with,  rotational  symmetry)  Into  which  a  supersonic 
stream  of  air  with  constant  velocity  may  enter  and  be  adla- 
batically  compressed  without  Interference  of  shocks*   In 
the  meantime  this  method  has  been  carried  out  in  a  simple 
case  and  the  results  are  now  presented  in  this  supplementary 
report. 

In  the  main  part  of  the  previous  report  a  method  was 
described  to  determine  isentroplc  flows  through  nozzles  with 
rotational  symmetry  by  employing  an  expansion  whose  terms  of 
order  zero  agree  with  the  results  of  the  hydraulic  treatment. 
Only  terms  of  first  and  second  order  were  then  considered 
for  the  refined  treatment.   To  check  the  accuracy  of  this  pro- 
cedure, terms  of  third  and  fourth  order  have  now  been  taken 
into  account  in  a  simple  case.   The  results  are  presented  in 
the  second  section  of  this  supplementary  report. 

Furthermore,  this  supplementary  report  contains,  in 
its  third  section,  a  more  detailed  derivation  of  the  basic 
differential  equations  employed. 
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1.   Construction  of  a  Supersonic  Compressor  or  Duct» 

For  the  construction  of  a  compressor  duct  It  has 
been  assumed  that  viscosity  and  heat  conduction  may  be 
neglected.   In  the  actual  design  of  a  compressor  duct 
this  assumption  can  no  longer  be  upheld  and  for  this 
reason  the  design  presented  here  can  be  considered  only 
an  approximation.   Nevertheless,  this  approximation  prob- 
ably gives  a  very  good  indication  of  the  proper  kind  of 
contour  shapes,  and  in  particular  of  the  length  of  the  con- 
verging section  of  the  duct  for  a  prescribed  compression  ratio. 

Since  viscosity  and  heat  conduction  are  ignored,  the 
flow  through  a  supersonic  cor^ressor  duct  is  the  reverse  of 
the  flow  through  a  perfect  nozzle  (a  nozzle  being  called 
perfect  if  it  produces  an  exhaust  flow  with  constant  velocity). 
For  convenience  we  shall  in  the  following  refer  to  exhaust 
flow  rather  than  to  compressor  flow. 

The  existence  of  perfect  three-dimensional  nozzle  flows 
had  been  indicated  by  A.  Busemann  [15],  1942^'"  who,  however, 
did  not  describe  a  construction.   It  had  escaped  the  attention 


*  Also  for  another  reason  the  presented  design  is  only  an 
approximation,  viz.,  because  shocks  must  be  admitted  in 
the  flow  since  otherwise  the  flow  would  be  unstable. (See 
Section  9  in  the  previous  report. ) 

-;:-»-  Reference  numbers  refer  to  the  bibliography  of  the  previous 
report. 
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of  the  writer  of  the  previous  report  that  a  construction 
of  a  perfect  nozzle  had  already  been  described  earlier 
by  P.  J.  Frankl  [21],  1934,  A  detailed  outline  of  the 
construction  of  a  "Prankl  nozzle"  (consisting  of  22  steps) 
is  described  by  M.  S.  Kisenko,  (see  bibliography  at  the 

end  of  this  report). 

Other  references  concerning  nozzle  flow  which  in 
the  meantime  came  to  the  attention  of  the  writer  are  listed 
in  the  bibliography. 

The  design  described  in  the  present  report  employs 
the  steps  outlined  in  Appendix  II,  Section  13  of  the  pre- 
vious report  with  minor  modifications  which  will  be  indi- 
cated. 

We  denote  the  pressure  by  p,  the  density  t>jf>,    the 
flow  speed  by  q,  the  sound  speed  by  c,  the  Mach  number  q/c 
by  M,  the  cross-section  area  by  A.   The  subscript  c  refers  to 
the  chamber,  -x-  to  the  throat,  e  to  the  mouth. 

A  pressure  ratio  Pq/Pq  =  5.66  was  assumed;"  it  corre- 
sponds to  the  expansion  ratio  A^/A^.  =  1.43  and  the  exhaust 
Mach  number  M^  =  1»79.   See  the  first  line  in  Table  I. 

A  regular  non -perfect  exhaust  flow  P^  was  constructed 

*by  the  method  of  Part  I  of  the  previous  report.   The  point 

A  on  the  axis  was  found  at  which  the  desired  Mach  number 
o 

M  was  reached.  A  Mach  line  0^  in  backward  direction  through 


-»-  A  value  somewhat  less  than  6  was  desired;   the  value  5.66 
was  convenient  for  numerical  computation. 
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A  was  determined.  (See  Section  7  of  the  previous  report.)  The 
flow  F  was  cut  off  beyond  the  transition  line  C  ,  and  there 
replaced  by  a  flow  G  which  is  turned  over  Into  a  parallel  flow 
(see  Sections  7  and  13).   The  results  are  shown  In  Table  II. 
This  figure  exhibits  firstly  the  streamlines  07  =  const.   The 

quantity  07  Is  the  distance  of  the  streamline  from  the  axis  at  the 

2 
throat*   The  quantity  Tr^^,q.„9^   Is  then  the  flvix  of  mass  per 

xinlt  time  through  the  tube  bounded  by  the  surface  7)  =   const. 

Further,  the  lines  q  =  const,  or  what  Is  equivalent, 
p  =  const,  or  M  =  const,  are  shown.   They  are  labeled  by  the 
values  of  the  Mach  niimber  M  along  these  lines.   The  correspond- 
ing values  of  P^/p  and  other  relevant  quantities  constant  on 
these  lines  are  given  in  Table  I.   The  line  M  =  1  separates 
the  region  of  subsonic  and  supersonic  flow.   The  line  M  =  1.79 
is  the  straight  Mach  line  D^  starting  at  the  point  A  on  the 
axis  beyond  which  the  flov/  has  constant  velocity  q  in  the 
axial  direction. 

Finally,  a  set  of  Mach  lines  is  drawn  in  the  supersonic 
range . 

One  observes  that  the  backward  Mach  lines  C  (and  also 
the  lines  M  =  const. )  intersect  each  other  near  the  throat 
cross-section  for  about  97  =  »7.   Continued,  these  lines  would 
have  a  cusp-like  envelope,  a  "limiting  line";   and  the  field 
of  flow  would  fold  over.   Clearly,  the  flow  could  not  be  used 
beyond  the  streamline  7^  =  .6  and  even  the  location  of  the  line 
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7^,  =  .6  is  not  quite  certain. 

Each  of  the  streamlines  shown  in  the  figijre  may 
serve  as  a  nozzle  contour.  For  the  purpose  of  comparison 
these  different  streamlines  were  so  magnified  that  they 
correspond  to  the  same  throat  diameter  and  then  superim- 
posed.  The  result  is  shown  in  Table  III,   The  coordinates 
of  points  on  these  contours  corresponding  to  throat, 
transition,  and  exit  are  given  in  Table  IV;  R  being  the 
radius  of  curvature  at  the  throat.   The  contour  correspond- 
ing to  77  =  ,1  is  omitted  because  it  would  be  too  long;   even 
the  contour  t^  =  .2  is  rather  long.   It  appears  reasonable  to 
employ  the  line  7^  =  .5  for  design.   (The  line  77  =  ,6  is  not 
much  shorter  and  its  location  is  uncertain. )   The  contoior 
97  =  .5  yields  a  length  of  2,5  throat  diameters  between  throat 
and  exit,  a  length  which  appears  rather  acceptable. 

As  to  the  details  of  the  calculation,  the  following 
may  be  mentioned.   The  basic  flow  F^  was  chosen  as  one  whose 
streamlines  are  essentially  hyperbolas;   more  exactly  h  =/l+|^ 
in  the  terminology  of  Part  I  of  the  previous  report.   The 
streamlines  were  dra;vn  according  to  formulas  (2. 11)  and  (2.12), 
the  latter  formula  being  simplified  to  y  =  h  77  since  the 
third  order  term  turned  out  to  be  very  small. 

The  Mach  line  C^  was  drawn  from  the  differential  equa- 
tion (see  13.09  and  13.10) 

dx  +  tan{/3   +  9)dy  =  0, 
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Table  III 

A  set  of  perfect  nozzle  contours  with 
the  same  throat  derived  by  magnifying 
the  streamlines  (97  =  constant)  of  the 
flow  shown  in  Table  II,   The  points  of 
inflection  are  indicated. 
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A  Bet  of  perfect  noztle   oontoura  with 
the  same   throat  derived  by  macnlfying 
the   streamiinoe    ( i?  -constant)  of  the 
flow  shown  in  Table   '1.     The  points  of 
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^   =  90°- a  here  being  the  angle  which  was  called'??  in 
Section  13,  and  oc  being  the  Mach  angle.  Expansions  of 
3C,  y,  ©»/3  with  respect  to  9p  were  inserted  and  terras  of 
second  and  higher  order  were  neglected.   The  result  was 
the  equation 

d^  +  li^j(^)tan/3^(  \  )d7^  =  0. 

Here  90°  -  f3     is  the  Mach  angle  corresponding  to  q/q.«.  =  Qq^^  ^» 
hence  cos  {5^   =  M  can  be  expressed  in  terms  of  qQ(  ^  )• 
The  equation  leads  to  the  integral 


d| 


^o^  ^  )tan/3o(  ^  J 

which  could  be  evaluated.  Since,  however,  just  the  points 
with  7^  =  .1,  ,2,  .*  were  desired  in  particular,  a  finite 
difference  procedtire  was  enrployed  with  the  step  A>7  =  •!• 
The  Mach  line  C  could  then  be  drawn. 

The  flow  denoted  by  G  in  Section  7  could  then  be 
determined  by  the  method  indicated  in  Section  13.   The  only 
modification  was  that  the  angle  yS  (there  called  99  )  as  a 
function  of  t,   was  not  read  off  from  Busemann's  diagram  but 
obtained  by  inversion  from  a  table  giving  £>  as  a  function 
of  yd  easily  obtained  on  the  basis  of  (13.03)  and  (13.02). 
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We  chose  the  parameters  <r  and  r   (see  Pig.  16  in  the 
previous  report)  so  that  they  were  identical  with  the  quantity 
ri    on  the  Mach  lines  D^  and  C^  respectively;  Ar  =  /^r  =  .1 
was  chosen  as  the  step.   The  advantage  of  doing  so  was  that 
then  the  intersections  of  the  streamlines  17  =  .  1,  .2,  ..  with 
these  Mach  lines  were  just  the  net  points  on  them. 

The  procedure  indicated  in  Section  13  then  proved  to 
he  manageable.   Near  the  axis  it  might  perhaps  have  been 
simplified  since  it  turned  out  that  the  Mach  lines  D  were 
nearly  straight  there.   The  directions  of  these  lines  on  the 
Mach  line  C  being  given  from  the  main  flow  P^,  one  might 
have  drawn  these  Mach  lines  straight  and  then  proceeded  in 
the  same  manner  as  for  the  corresponding  two-dimensional  con- 
struction. 

2.   The  Accuracy  of  the  Determination  of  Nozzle  Flow  by  the 
Method  of  Expansion. 

To  obtain  some  information  on  the  accuracy  of  the  method 
for  determining  the  nozzle  flow  explained  in  Part  I  of  the 
previous  report,  the  expansion  has  now  been  carried  through 
up  to  the  fourth  order  in  a  simple  example  and  the  results 
have  been  compared  with  those  obtained  if  only  terms  up  to 

the  second  order  are  considered. 

The  example  taken  is  the  same  as  that  used  for  the  basic 

flow  P  in  the  construction  of  the  preceding  section,  viz.: 
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yleldlng  streamlines  which  are  approximately  hyperbolas. 
The  results  of  drawing  streamlines  from 

are  compared  with  those  resulting  from 


^  =  i  ■*■  ^2'^^        '   7  =  M 


in  Table  V  .   It  is  seen  that  deviations  are  noticeable  for 
1^  =  .5  and  77  =  •&,      Up  to  7^  =  .4  the  first  approximation 
might  be  considered  sufficient. 

The  coefficients  Xg,  yg  are  given  in  (2.11)  and  (2.12). 
The  coefficient  x^  is 


with 


Qo  =  M^   - 

^O              0 

2  ''o 

1   , 

^1    -  J^' 

For  hQ(  ^  )  =  /l  +  ^^  »  ir^  particular 


1  _j „2    K/q^)^  -^  ^ 

^4     16     V  2   o   ,   ,   ,2 
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3»  Derivation  of  the  Formulas  In  Appendix  I  of  the  Previous 
Report. 

It  had  been  suggested  that  derivations  be  given  for 
the  formulas  (12.01)  In  Appendix  I.   These  formulas  give 
expressions  for  the  derivatives  of  x  and  7  with  respect  to 
I  and  7^  ,      Of  course,  they  can  be  derived  formally  from 
the  equations  (2.01),  (2.02)  which  give  expressions  for  the 
derivatives  of  fi   and  with  respect  to  x  and  y,  observing  that 
^  Is  a  function  of  ^  by  (2.09)  and -j^  a  function  of  7^  by 
(2.10).   In  the  following  equation  (12.01)  Is  derived  directly 
from  the  continuity  relation  and  the  condition  that  the  flow 
be  Irrotatlonal. 

The  two  partial  derivatives  ~   ,  ^  Imply  that  In  carry- 
ing out  the  differentiation,  the  quantity  yr  should  be  kept 
constant.   Since  the  equation  y^  =  const.  Is  that  of  a  stream- 
line. It  Is  clear  that  -^   ,  -^  are  derivatives  along  a  stream- 
line.  Consequently  these  two  quantities  represent  the  component 
of  a  vector  which  Is  tangential  to  the  streamline  through  the 
point  In  question.   In  other  words,  this  vector  points  In  the 
same  direction  as  the  velocity  vector.   The  angle  the  velocity 
vector  makes  with  the  axis  was  denoted  by  9.   Let  further  s  be 
the  arc  length  along  this  streamline.   Then  we  have 


^  =  d3     Q    ^  =  ds 
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where 

ds 
djif 


M  *{W- 


ds 
The  factor  -r—  can  easily  be  evaluated.   Prom  the 

definition  of  the  potential  function  0  we  have 

-|  =  q  =  |/u  +  V   , 

hence  ^  =  1/q  and 

^  -   QO°  Q  9y  _     sin  Q 

ajrf  q  •      ajjJ  q        • 

'iflfhen  i   is   introduced  as  new  variable   instead  of  jtf  by 
d0  =  q^d^     ,    these  relations  become 

__  =  -ii.  cos  0  ,     -Tt-  =  —  sin  9     • 

The  partial  derivatives  -^  ,  ^  are  correspondingly 
the  con5)onents  of  a  vector  tangential  to  the  potential  line 
through  the  point  in  question.  Hence 

^  =  -  -^  sin  9  ^  =  ^  cos9  . 

ay    dV' ^^"^  ^'  av^  dy' '^°^®  • 

Here   s   is   the  arc   length  along  the  potential  line,   and 


ds   ^    //ax\2  /^^^  2 

dV     Y    XdYl      ""    laW 
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We  assume  that  s  Increases  away  from  the  axis;   then  on  the 

axis  dy  =  ds;   therefore  the  signs  are  correct  as  indicated. 

ds 
To  evaluate  -^  we  note  that  according  to  the  definition 

of  the  stream  function 

d{vy^)   =    yoqd(2Trs)  , 

the  left  hand  side  being  the  mass  fliuc  per  unit  time  crossing 
between  two  stream  tubes  y =  const,  for  which  the  difference 
in  the  values  of  y-  is  dyr •      The  right  hand  side  is  the  product 
of  the  mass  flux  per  unit  time   per  unit  area  c  q  and  the 
area  of  the  section  of  the  potential  surface  cut  out  by  these 
t?/o  stream  tubes.  Prom  this  relation  one  has 

•^dy'=  ^  qds 


or 


Hence 


da  _  yr 


^=  -  J:  Sine,  |l  =  XcosG 
^Y  yoq        9r  ^q 


When  TO  is  introduced  as  new  variable  by 
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then  these  equations  become 


|S  =   _  h^sin  0  ,   |2:  =  h^cos  9  , 


where  the  abbreviation 


2   /°^fr^-::- 
h^  =  '^ 


/oq 


has  been  used. 
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